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Abstract 

Three  dimensional  stress  analysis  of  spiral  bevel  gears  in 
mesh  using  the  finite  element  method  is  presented.  A 
finite  element  model  is  generated  by  solving  equations  that 
identify  tooth  surface  coordinates.  Contact  is  simulated 
by  the  automatic  generation  of  non-penetration  constraints. 
This  method  is  compared  to  a  finite  element  contact 
analysis  conducted  with  gap  elements. 

Introduction 


The  research  reported  in  this  study  will  utilize  the 
numerical  solution  for  spiral  bevel  surface  geometry  with 
the  automatic  generation  of  non-penetration  constraints  to 
model  the  contacts  during  meshing.  Pinion  tooth  and  gear 
tooth  surfaces  will  be  developed  based  on  the 
manufacturing  kinematics.  A  multi-tooth  model  (4  gear 
and  3  pinion  teeth)  is  created  by  replicating  pinion  and 
gear  teeth  in  space  with  a  solid  modeler  [10].  The  results 
for  the  gap  element  model  are  compared  to  the  non¬ 
penetration  constraint  model. 


Spiral  bevel  gears  are  used  to  transmit  power  between 
intersecting  shafts.  In  helicopter  transmissions  these  gears 
operate  at  very  high  speeds  and  transmit  substantial 
power.  To  reduce  vibration,  noise  and  in^jrove  life,  the 
stresses  and  deflections  during  meshing  must  be 
understood. 

The  finite  element  method  (FEM)  is  particularly  well 
suited  to  study  tooth  deflections,  contact  stresses  and  tooth 
fillet  stresses.  Much  effort  has  focused  on  predicting 
stresses  in  gears  with  FEM.  Most  of  this  research,  to 
date,  has  been  done  on  parallel  axis  gears  with  two 
dimensional  geometry.  Because  of  the  difficulty  of 
correctly  identifying  the  three  dimensional  geometry,  only 
a  few  researchers  have  investigated  spiral  bevel  gears  with 
FEM  [1,2]. 

Previously  a  single  spiral  bevel  tooth  has  been  analyzed 
with  FEM  [3],  however  in  that  analysis  a  contact  stress 
distribution  was  assumed.  To  model  the  correct  contact 
load  distribution  requires  a  gear  and  pinion  loaded  in 
mesh.  In  reference  [4],  a  gearset  in  mesh  was  analyzed 
using  gap  elements  to  model  the  contact  zones. 


Tooth  Surface  Coordinate  Equations 

The  system  of  equations,  required  to  dehne  the  coordi¬ 
nates  of  spiral  bevel  gear  tooth  surfaces,  were  derived  in 
reference  3,  and  are  briefly  summarized  here. 

The  first  equation  (equation  of  meshing)  is  based  on  the 
kinematics  of  manufacture  and  the  machine  tool  settings. 
The  equation  of  meshing  requires  the  relative  velocity 
between  a  point  on  the  cutting  surface  and  the  same  point 
on  the  pinion  being  cut  must  be  perpendicular  to  the 
cutting  surface  normal. 

n*  V  =  0 

where, 

n  is  the  normal  vector  from  the  cutter  surface  and 
V  is  the  relative  velocity  betweoi  the  cutter  and  the 
workpiece  surfaces  at  the  specified  location.  This 
equation  is  developed  in  terms  of  machine  tool  parameters 
U,  0  and  where,  U  and  0  (length  and  rotational 
orientations)  locate  a  point  on  the  cutting  head  cone  and 
is  the  rotated  orientation  of  the  cutter  as  it  swings  on 
the  cradle. 
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Since  the  kinematic  motion  of  cutting  a  gear  is  equivalent 
to  the  cutting  head  meshing  with  a  simulated  crown  gear, 
an  equation  of  meshing  can  be  written  in  terms  of  a  point 
on  the  cutting  head  i.e,  in  terms  of  U,  6  and  The 
equation  of  meshing  for  straight  sided  cutters  with  a 
constant  ratio  of  roll  between  the  cutter  and  the  workpiece 
is  given  by  (ref.  5): 

(U-r  cotijr  cosi|r)  cosy  sint 
+  S(m^  -  siny)  cosijr  sm0 
T  cosy  siny  sin(q  -  <t»^  (1) 

±  (cosy  siny  +  siny  cosy  cost) 

-  siny  cosi|r  sinT  =  0 

The  upper  and  lower  signs  are  for  left  and  right  hand 
gears  respectively.  The  following  machine  tool  settings 
are  defined: 

Ip  cutting  blade  angle 

T  (0Tq±^j) 

q  cradle  angle 

7  root  angle  of  workpiece 
machining  offset 

L„,  vector  sum  of  change  of  machine  center  to  back 
and  the  sliding  base 

“cw  *11®  relationship  between  the  cradle  and  the 

workpiece  for  a  constant  ratio  of  roll 

U  generating  cone  surface  coordinate 

S  radial  location  of  cutting  head  in  coordinate 
system  S„ 

r  radius  of  generating  cone  surface 

Equation  (1)  is  equivalent  to: 

/i  iU,  0,  ♦;  =  0  (2) 

Because  there  are  three  unknowns  U,  0,  and  three 
equations  must  be  developed  to  solve  for  the  surface 
coordinates  of  a  spiral  bevel  gear.  The  three  parameters 
U,  0  and  are  defined  relative  to  the  cutting  head  and 
cradle  coordinate  systems  (Sc  and  SJ  respectively.  These 
parameters  can  be  transformed  through  a  series  of 
coordinate  transformations  to  a  coordinate  system  attached 
to  the  workpiece  [5].  Or  U,  0  and  ip^  can  be  mapped 
into  X*,  Y»,  in  coordinate  system  S„  attached  to  the 
workpiece.  These  transformations,  used  in  conjunction 
with  two  other  geometric  requirements,  give  the  two 
additional  equations. 

The  correct  U,  0  and  4>c  *li2*  solves  the  equation  of 
meshing,  must  also,  upon  transformation  to  the  workpiece 
coordinate  system  Sw>  result  in  a  axial  coordinate  that 
matches  with  the  value  of  Z  found  by  projecting  the  tooth 


face  in  the  XY  plane. 

-  Z  =  0  (3) 

This  equation  along  with  the  correct  coordinate  trans¬ 
formations  (see  Equation  10)  result  in  a  second  equation 
of  the  form: 

f2(u,  e,  ♦c )  =  0 

A  similar  requirement  for  the  radial  location  of  a  point  on 
the  workpiece  results  in  the  following: 

r  -  ^(Xl  +  =  0 

These  projections  are  as  shown  in  figure  1.  The 
appropriate  coordinate  transformations  (see  Equation  10) 
will  convert  equation  6  into  a  function  of  U,  0  and  <p^ 

MU,  0,  ♦;  =  0  (6) 

Equations  (2),  (4),  (6)  form  a  system  of  nonlinear 
equations  necessary  to  define  a  point  on  the  tooth  surface. 

Solution  Technique 

Equations  (2),  (4),  (6)  form  a  nonlinear  system  of  equat¬ 
ions  that  do  not  have  a  closed  form  solution.  They  are 
solved  using  Newton’s  method  (ref.  8). 

An  initial  guess  Ug,  0»,  ^co  >s  used  to  the  start  iterative 
solution  procedure.  Newton’s  method  is  used  to 
determine  subsequent  values  of  the  updated  vector  (U,,, 

0k»  0ck)- 


U, 

ry’ 

^t-1 

e* 

= 

0*-i 

+ 

rU 

kJ 

where  the  vector  Y  is  the  solution  of: 

The  3x3  noatrix  in  the  preceding  equation  is  the  Jacobian 
matrix  and  must  be  inverted  each  iteration  to  solve  for  the 


^i(e*-’)  ^i(<t.r') 

dU  00  04)^ 

if'* 

fi(l7‘-«.0*-‘,4>f') 

^2(t/*-‘)  ^2(e*-‘)  ^2(<i>r') 

•if'* 

4(l/‘-*.0‘-*,4»f*) 

dU  00  04), 

^3(1/*-*)  ^3(0*-*)  ^3(4>*’‘) 

dU  00  04), 

if'* 

f3(l/‘-*.0*-*,4>^‘) 

(8) 
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Y  vector.  The  equation  of  meshing,  function  fj,  is 
numerically  differentiated  directly  to  find  the  terms  for  the 
Jacobian  matrix.  Function  f2  and  f3  cannot  be  directly 
differratiated  with  respect  to  U,  0  and  After  each 
iteration  0“,  (“  cutting  head  coordinate 

system)  are  transformed  into  the  woiicpiece  coordinate 
^stem,  S„,  with  the  series  of  coordinate  transformations 
as  givra  in  Equations  9  and  10.  (The  right  side  of 
equation  9  locates  a  point  on  the  cutting  head  cone  as 
shown  in  figure  2.) 


K 

r  cotijr  -  V  cosijr 

=  [3fJ 

V  sini|r  sin6 

1]  simlr  cos6 

1 

1 

where, 

[Af^] =[M^/(<|>^][Af (10) 

The  matrix  [M*J  above  represents  a  series  of 
homogeneous  coordinate  transformations  from  a 
coordinate  on  the  cutting  head  onto  a  coordinate  system 
on  the  workpiece  [5]. 

Functions  f,  and  fj  are  evaluated  by  starting  with  an  initial 
Ufc,  Ok  and  performing  the  transformations  in 
Equation  (10)  and  evaluating  Equations  (4)  and  (6).  The 
numerical  differentiation  of  and  f3  is  performed  by 
transforming  Uk+inc,  Ok+inc,  ^k'l'inc  (where  inc  is  a 
small  increment  iqrpropriate  for  numerical  differentiation) 
into  X,+inc,  Y,+inc,  Z*+inc.  Equations  (4)  and  (6) 
are  thra  used  to  evaluate  the  numerical  differentiation. 
Function  f^,  fj,  f3  and  their  partial  derivatives  are  required 
for  the  Jacobian  matrix  and  are  updated  each  iteration. 
The  iteration  procedure  continues  until  the  Y  vector  is  less 
then  a  predetermined  tolerance.  This  completes  the 
solution  technique  for  a  single  point  on  the  spiral  bevel 
gear  tooth  surface.  In  this  way  the  coordinates  of  the 
surface  of  the  tooth  are  defined.  This  solution  technique 
is  repeated  for  each  of  the  four  surfaces;  gear  convex, 
gear  concave,  pinion  convex  and  pinion  concave. 

Since  all  four  sur&ces  are  generated  independently, 
additional  matrix  transformations  are  required  to  obtain 
the  correct  tooth  thickness.  The  concave  surfaces  are 
fixed  on  each  tooth.  The  convex  surfaces  are  rotated  as 
required.  The  angle  of  rotation  is  obtained  by  matching 
the  tooth  top  land  thickness  with  the  desired  value. 


Gear  and  Pinion  Orientations  Required  for  Meshing 

When  generating  the  pinion  and  gear  surface  as  described 
above,  the  pinion  cone  and  gear  cone  apex  will  meet  at 
the  same  point.  This  point  is  the  origin  of  the  fixed 
coordinate  system  attached  to  the  workpiece  being 
generated.  To  place  the  gear  and  pinion  in  mesh  with 
each  other,  rotations  described  in  the  following  exanqile 
are  required  (ref.  6). 

1.  The  gear  tooth  surface  points  are  rotated  by 
360/N,+ 180  CW  about  the  global  2^  axis,  for  this 
example,  the  rotation  is  190  degrees. 

2.  The  pinion  is  rotated  by  90  deg  CCW  about  the 
global  Y  axis. 

3.  Because  the  four  surfaces  are  defined 
independently,  their  orientation  is  random  with 
respect  to  meshing.  The  physical  location  of  the 
gear  and  pinion  after  rotations  described  above 
correspond  to  the  gear  and  pinion  overlapping. 
To  correct  this  condition  the  pinion  is  rotated  CW 
about  its  axis  of  rotation  xmtil  surface  contact 
occurs.  For  the  example  used  in  this  study,  the 
rotation  was  3.56  deg. 

To  make  a  complete  pinion,  the  generated  pinion  tooth 
was  copied  and  rotated  12  times  and  the  generated  gear 
tooth  was  copied  and  rotated  36  times  with  the  aid  of  a 
solid  modelling  program  as  shown  in  figure  3. 

Finite  Element  Analysis  of  Spiral  Bevel  Gears 

Recent  finite  element  analysis  research  on  spiral  bevel 
gears  has  attempted  to  solve  the  contact  stress  distribution 
in  a  multi-tooth  model  (4  gear  and  3  pinion  teeth)  (ref.  4). 
The  tooth  pair  contact  zones  in  reference  3  were  modeled 
with  gap  elements.  In  the  current  study,  use  of  software 
with  automated  contact  options  is  used.  This  is  done  to 
avoid  certain  limitations  in  the  use  of  gap  elements,  such 
as: 

(i)  Difficulties  in  connecting  the  gap  elements  with 
proper  orientations. 

(ii)  The  orientation  of  the  contact  plane  remains 
unchanged  during  deflection. 

(iii)  Difficult  to  accurately  select  the  properties  such 
as  appropriate  open/closed  stiffness  values. 
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selection  of  the  stifhiess  matrix  update  strategy 
and  efficient  problem  restarts. 

New  advances  in  the  state  of  the  art  for  FEM  provide 
deformable  body  against  deformable  body  penetration 
algorithms  udiich  can  be  used  to  establish  the  nonlinear 
boundary  conditions  for  contact  problems.  MARC  (ref. 
9)  is  one  such  FEM  package  software  which  uses  this 
algorithm  to  automatically  detect  contact. 

Automated  Contact  Analysis 

Situations  in  which  contact  occurs  are  common  to  many 
different  nonlinear  applications.  Contact  forms  a 
distinctive  and  important  subset  to  the  category  of 
changing-status  nonlinearities.  Contact,  by  its  very  nature, 
is  a  nonlinear  problem.  During  contact,  both  the  forces 
transmitted  across  the  area  and  the  area  of  contact  change. 
The  force-displacement  relationship  thus  become 
nonlinear.  Usually,  the  transmitted  load  is  in  the  normal 
direction. 

The  program  used  MARC  (ref.  9),  is  a  general  purpose 
computer  program  for  linear  and  nonlinear  stress  and  heat 
transfer  analysis.  This  program  is  capable  of  solving 
problems  with  nonlinearities  that  occur  due  to  material 
properties,  large  deformations,  or  boundary  conditions. 
In  general,  the  solution  of  nonlinear  FEM  problems 
requires  incremental  solution  schemes  and  sometimes 
requires  iterations  within  each  load/time  increment  to 
satisfy  equilibrium  conditions  at  the  end  of  each  step. 

The  FEM  program  used  has  a  fully  automatic  contact 
option  which  enables  the  analysis  between  finite  element 
bodies  without  the  use  of  any  special  gap  or  contact 
elements.  The  contact  option  was  originally  designed  for 
aiudysis  of  manufacturing  processes  such  as  forging  or 
sheet  metal  forming,  but  its  capabilities  have  been 
expanded  to  meet  other  analysis  requirements.  The  work 
presented  here  utilized  the  program  of  reference  9  running 
on  a  supercomputer. 

Contact  between  the  bodies  is  handled  by  imposing  non¬ 
penetration  constraints.  The  non-penetration  constraint, 
as  shown  in  figure  4  and  figure  S  is: 

Ux  •  n  <  D 

where:  is  the  displacement  of  node  A, 

n  is  normal  to  the  contacted  surface, 

D  is  the  contact  tolerance 


Solver  constraints  are  used  to  impose  the  non-penetration 
constraint,  and  a  very  efficient  surface  contact  algorithm 
which  allows  the  user  to  simulate  general  2D/3D  multi¬ 
body  contact.  Both  "deformable-to-rigid"  and 
’deformable-to-deformable"  contact  situations  are 
allowed.  The  user  needs  to  only  identify  which  bodies 
might  come  in  contact  during  the  analysis.  The  bodies 
can  be  either  rigid  or  deformable,  and  the  algorithm 
tracks  variable  contact  conditions  automatically.  Thus, 
the  user  no  longer  needs  to  worry  about  the  location  and 
open/close  status  checks  of  ’gap  elements*.  Automatic, 
in  this  context,  means  that  user  interaction  is  not  required 
in  treating  multibody  contact  and  friction,  and  the 
program  has  automated  the  imposition  of  non-penetration 
constraints. 

Real  world  contact  problems  between  rigid  and/or 
deformable  bodies  are  actually  3D  in  nature.  To  solve 
such  contact  problems,  one  needs  to  define  bodies  and 
their  boimdary  surfaces. 

Deformable  bodies  are  defined  by  the  elements  of  which 
they  are  made.  Once  all  the  boundary  nodes  for  a 
deformable  body  are  determined,  4  point  patches  are 
automatically  created,  which  are  constantly  updated  with 
the  body  deformation.  Contact  is  determined  between  a 
node  and  all  body  profiles,  deformable  or  rigid. 

The  user  must  define  bodies  so  that  their  boimdary  sur¬ 
faces  can  be  established.  Deformable  bodies  are  defined 
by  a  list  of  finite  elements.  Because  the  contact  boundary 
conditions  are  a  ftmction  of  the  applied  load,  the  analysis 
must  be  carried  out  incrementally.  Within  each  load  or 
time  increment  of  an  analysis,  additional  iterations  may  be 
required  to  stabilize  the  contact  zone.  Contact  problems 
involve  two  important  aspects: 

(i)  the  opening  and  closing  of  the  gap  between  bodies 

(ii)  friction  between  the  contacting  surfaces. 

The  MARC  program  establishes  a  hierarchy  between  the 
bodies  so  that  at  a  given  contact  interface,  one  body  is  the 
contactor  and  the  other  body  is  the  contacted.  The  set  of 
nodes  on  the  boundary  surface  of  a  contactor  are 
candidate  nodes  for  contact.  The  boundary  surface  of  a 
contacted  body  is  defined  by  a  set  of  geometrical  entities. 
A  user  specified  contact  tolerance  is  used  to  determine  the 
body  separation  distance  which  determines  whether  two 
bodies  are  considered  to  be  in  contact  with  each  other. 
The  contactor’s  boundary  nodes  are  prevented  from 
penetrating  the  surface  of  the  contacted  body  by  imposing 
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solver  constraints.  For  contact  between  two  deformable 
bodies,  MARC  applies  multipoint  constraints  in  the  form 
of  ties.  The  ties  link  the  motion  of  one  node  on  the 
contactor  body  to  two  adjacent  nodes  on  the  contacted 
body.  During  each  iteration  as  nodes  enter  and  leave 
contact  or  slide  betwe^  adjacent  node  segments  on 
contacted  bodies,  a  bandwidth  optimization  is  performed 
to  reduce  the  computer  processing  time  required. 

A  static  analysis  of  two  bodies  that  are  not  initially 
connected  poses  special  problems  with  a  FEM,  if  one  of 
the  bodies  has  a  rigid  body  motion  component.  If,  at  any 
time,  the  two  bodies  are  disconnected  then  the  sd^ess 
matrix  would  become  singular  and  unsolvable  (in  a  static 
analyses).  In  order  to  overcome  this  difficulty,  the  two 
gears  are  connected  with  weak  springs.  The  spring 
sdf&ess  must  be  negligible  compared  to  the  material 
sdfhiess. 

Model  descriptions 

Pinion  and  gear  design  data  and  generation  machine 
settings  are  as  shown  in  Tables  I  and  n.  This  information 
is  used  as  input  into  a  computer  program  that  solves  the 
3  nonlinear  equadons  (previously  described)  used  to 
define  the  points  on  the  pinion  and  gear  surface.  The 
program  goierates  a  NxM  mesh  on  the  tooth  profile  of  a 
^iral  bevel  gearset  [7].  After  one  gear  and  one  pinion 
tooth  was  made,  the  geometric  modeler  [10]  was  used  to 
rotate  the  teeth  to  create  4  gear  teeth  and  3  pinion  teeth  in 
mesh.  Eight  noded,  isoparametric  3D  brick  elements 
were  used.  The  seven  tooth  model  is  as  shown  in  figure 
6.  The  seven  tooth  model  consisted  of  8793  elements  and 
11261  nodes. 

Loading  and  Boundary  Conditions 

The  boundary  condidons  for  the  seven  tooth  model  are 
also  shown  in  figure  6.  The  boundary  condidons  are 
applied  such  that  the  gear  teeth  are  made  to  pivot  about  a 
fixed  point  along  the  axis  of  rotation  for  the  gear  on  the 
Z  axis.  The  nodes  where  the  forces  are  applied  are  in  the 
hub  of  the  gear.  Fixed  displacement  boundary  conditions 
are  applied  at  8  nodal  points  in  the  Z  direction  only  for 
the  gear  and  in  all  direcdons  for  the  pinion  as  shown  in 
figure  6.  Since  this  is  a  static  problem  involving  two 
bodies  (the  pinion  and  the  gear)  in  contact,  as  described 
in  a  previous  chq>ter,  weak  springs  are  added  to  prevent 
die  rotadonal  rigid  body  modes  for  the  gear.  Eight 
springs  are  added  away  from  the  contact  region.  The 
wrings  connect  the  comer  nodes  of  the  pinion  with  the 
comer  nodes  of  the  gear  on  the  faces  where  contact 


occurs.  The  sdffhess  of  the  springs  are  100  Ibs/in. 
This  is  insignificant  when  compared  to  the  tooth  contact 
sdf&ess  and  therefore  does  not  effect  the  overall  soludon. 
The  maximum  torque  for  the  gear  mesh  s&died  was  9508 
m.*lbs.  on  the  gear.  This  torque  is  applied  as  a 
concentrated  force  with  a  moment  arm  on  the  gear  hub. 
This  concentrated  force  for  the  seven  tooth  model  was 
4714  lbs.  with  a  moment  arm  of  2.017  mches.  The 
force  is  applied  mcrementally  for  convergence  to  occur. 

Gap  Element  and  Automated 

Contact  Analysis  Comparison 

Contact  stresses  on  spiral  bevel  gears  were  studied  by 
researchers  with  gap  elements  m  references  11  and  12  on 
a  similar  seven  tooth  model,  comparison  of  the  results 
with  automated  contact  analysis  is  presented.  Both 
models  contamed  the  same  mesh  density,  boundary 
conditions,  material  properties  and  loadmg.  The  nodal 
stress  results  of  pinion  tooth  #1  obtamed  from  automated 
contact  analysis  and  gap  elements  are  as  shown  m  figures 
7  and  8,  respectively.  &  both  analyses  the  highest 
concentrated  stress  value  occur  at  the  same  node. 
Comparison  of  the  results  of  these  two  runs  are  as 
follows. 


Gap  Element 

Automated 
Contact  Analysis 

Min  nodal  principal 
stress 

-296,410  (jpsi) 

-291,503  (psD 

CPU  time  (approx) 

30  min. 

80  min. 

Elemental  principal 
stress 

-84,761  (psi) 

-113,577  (psi) 

Gap  element  closed  or 
nodal  points  with  contact 

4 

8 

No.  of  iterations 

4 

8 

The  nodal  stresses  were  m  very  good  agreement  for  both 
the  gap  element  FEM  and  the  automated  contact  FEM 
model.  Elemental  and  nodal  stresses  differed  greatly  for 
both  models  mdicatmg  a  need  for  mesh  refinement.  The 
gap  element  model,  as  expected,  used  less  cpu  time. 
More  mformation  is  given  by  virtue  of  locatmg  and 
definmg  the  gap  elements.  This  reduced  the 
computational  time  required  to  search  for  and  calculate 
the  contact  area  and  stresses.  However,  as  noted  earlier, 
it  is  much  more  difficult  to  calculate  the  proper 
orientation  and  location  of  the  gap  elements.  It  is  more 
convenient  to  analyze  contact  with  the  automatic  contact 
FEM  program. 
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The  number  of  contact  nodes  at  the  contact  region  in  the 
automated  contact  FEM  program  was  higher  than 
identified  by  the  gap  element  FEM  program.  With  the 
pinion  considered  body  1  (contactor  body)  and  the  gear 
considered  body  2  (contacted  body),  eight  nodes  contacted 
as  shown  above.  With  body  1  and  body  2  switched,  16 
nodes  were  foimd  to  have  contact.  Presumably  this  sort 
of  discrepancy  occurred  because  the  mesh  was  too  coarse. 

CONCLUSIONS 

Three-dimensional  contact  analysis  of  spiral  bevel  gears 
in  mesh  was  performed.  Four  gear  teeth  and  three  pinion 
teeth,  generated  by  solving  equations  based  on  gear 
manufacturing  kinematics,  were  used  to  model  contact. 
An  automatic  contact  analysis  algorithm  is  utilized  with 
the  finite  element  method.  The  contact  algorithm  searches 
for  and  finds  contact  without  using  gap  elements. 

Comparison  of  nodal  and  elemental  principal  stresses  is 
favorable  with  a  previous  FEM  that  utilized  gap  elements. 
However,  for  both  models  the  nodal  and  elemental 
stresses  were  significantly  different.  Also,  both  models 
had  large  stress  gradients  at  the  contact  point.  This 
indicates  that  a  more  reEned  finite  element  mesh  is 
required. 

The  contact  algorithm  model  takes  more  cpu  time 
compared  to  the  gap  element  model.  The  contact 
algorithm  has  to  find  where  contact  occurs,  with  gap 
elements  contact  is  defined  to  occur  at  the  gap  elements. 

Although  the  gap  element  modal  ran  with  less  cpu  time, 
use  of  the  contact  algorithm  is  more  convenient.  The 
contact  algorithm  is  easier  because  gap  elements  do  not 
have  to  be  mathematically  defined  in  space  with  proper 
orientation. 
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TABLE  I:  PINION  AND  GEAR  DESIGN  DATA 


1  PINION 

GEAR 

Number  of  teeth  pinion 

12 

36 

Dedundum  angle,  degree 

1.5666 

3.8833 

Addendum  angle,  degree 

3.8833 

1.5666 

Pitch  angle,  degree 

18.4333 

71.566 

Shaft  angle,  degree 

90.0 

90.0 

Mean  spiral  angle,  degree 

35.0 

35.0 

Face  width,  mm  (in) 

25.4  (1.0) 

25.4  (1.0) 

Mean  cone  distance,  mm  (in) 

81.05  (3.19) 

81.05  (3.19) 

Inside  radius  of  blank,  mm  (in) 

5.3  (0.6094) 

3.0  (.3449) 

Top  land  thickness,  mm  (in) 

2.032  (0.080) 

2.489  (0.098) 

Clearance,  mm  (in) 

0.762  (0.030) 

0.929  (0.0366) 

TABLE  H:  GENERATION  MACHINE  SETTINGS 


When  contnet  is  detected,  Muhi-Puint 
Displacement  Constraint  Ties  are 
introduced  to  prevent  body  penetration 


Figure  5.  Contact  depiction 


Pinion 


Figure  6.  Seven  tooth  model  with  boundary  conditions 
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Figure 


A 


Figure  8. 


PSI 

A  85251. 
B  22459. 
C  -40334. 
D  -103126. 
E  -165918. 
F  -228711. 
G  -291503. 


7.  Nodal  stress  result  on  pinion  obtained  from 

automated  contact  analysis  in  seven  tooth  model 


PSI 


-17629. 

-57455. 

-97281. 

-137107. 

-176933. 

-216759. 

-256584. 

-296410. 


Nodal  Stress  Result  on  Pinion  Obtained  From  the  Gap 
Element  Solution  and  Seven  Tooth  Model 
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